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Introduction 


La mécanique des fluides est  l’étude du comportement 
des fluides (liquides et gaz) et des forces internes associées. C’est une branche 
de la mécanique des milieux continus qui modélise la matière à l’aide de 
particules assez petites pour relever de l’analyse mathématique mais assez 
grandes par rapport aux molécules pour être décrites par des fonctions 
continues. Elle se divise en deux parties : 

v la statique des fluides est l’étude des fluides au repos, qui se réduit pour 
l’essentiel à l’hydrostatique. 

v la dynamique des fluides, l’étude des fluides en mouvement 
(hydrodynamique). 

On distingue également d’autres branches liées à la mécanique des fluides : 

l'hydraulique, l'hydrodynamique, l'aérodynamique, ... Une nouvelle approche a 

vu le jour depuis quelques décennies: la mécanique des fluides numérique ou 

en anglais " Computational fluid dynamics". qui simule l'écoulement des 

fluides en résolvant les équations qui les régissent à l'aide d'ordinateurs très 

puissants : les supercalculateurs. 

L’étude de la mécanique des fluides remonte au moins à l’époque de la Grèce 

antique avec Archimède qui fut à l’origine de la statique des fluides. 

Dans certains problèmes particuliers, faute de modélisation numérique correcte 

des phénomènes, des modèles réduits sont utilisés. Pour cette raison, et aussi 

pour présenter des lois empiriques, la mécanique des fluides utilise 


systématiquement des nombres sans dimension. 


Ceci constitue le document de cours de Dynamique des fluides réels destiné 
aux étudiants de première année Master de l’université Dr Moulay Tahar de 
Saida ayant choisi le département de Génie des procédés. Nous limiterons notre 
étude à celle des fluides réels. Les lois et modèles simplifiés seront utilisés 


pour des fluides continus dans une description macroscopique. 


Ce cours est subdivisé en quatre chapitres : 

Le chapitre I est focalisé sur les aspects physiques du phénomène de viscosité et 
les lois de comportement des fluides (Newtonien, non newtonien,....). 

Le chapitre II est consacré à l’étude de la dynamique des fluides visqueux 
incompressibles dans laquelle on classe les écoulements laminaires et 
écoulements turbulents. La résolution de l'équation fondamentale de Navier- 
Stokes en dynamique des fluides passe par l'introduction de la notion de pertes de 
charge. Ainsi, pour rendre compte de la dissipation d'énergie due aux frottements 
visqueux, ces pertes de charges prendront place dans la formulation d'une équation 
de Bernoulli généralisée. 

Par rapport aux écoulements externes, les écoulements internes (écoulement 
parallèle dans un conduit bidimensionnel, écoulement de Couette, écoulement 
de Poiseuille et écoulement dans un espace annulaire) ont ceci de particulier 
que les conditions aux limites y sont omniprésentes. Ils sont abordés dans le 
chapitre III sous deux aspects complémentaires : les profils de vitesse et les 


pertes de charge qui conditionnent le dimensionnement des circuits fluides. 


A lissu du chapitre IV, les équations adimensionnelles et les conditions de 
similitude sont y énoncées et explicitées à partir de l'analyse dimensionnelle 


d'une relation mathématique entre grandeurs physiques. 
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I-1Mise en évidence : viscosimètre de Couette 

Le viscosimètre (Figure I-1)est constitué de deux cylindres coaxiaux l'un est mis en 
rotation avec la vitesse angulaire œ. L'autre est entraîné par les forces de viscosité et est 
maintenu immobile par l'application d'un couple C. L'espace entre les deux cylindres est 
faible et rempli du fluide dont on veut mesurer la viscosité. 

On appelle R; et R, les rayons des cylindres et e la différence R — Rize. e est très 
petite de telle façon que l'on puisse considérer la répartition des vitesses comme linéaire 
entre les deux cylindres. 


On montre que la force de frottement est proportionnelle au gradient de vitesse soit : 


SV 
F = p70) 


Cylindre maintenu 
immobile 


Fluide 


Cylindre en rotation 





Figure I-1 : Viscosimètre de Couette. 


Pour définir une grandeur capable de quantifier la notion de viscosité, on recherche une 
relation mathématique entre la force de déformation et la vitesse du fluide Soit une 


couche de fluide placée entre deux plaques (Figure I-2) 


Plaque mobile 





Figure I-2 :Couche de fluide placée entre deux plaques. 
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Remarques 

v Lorsque la vitesse n’est pas trop grande, les couches de fluide ne se mélangent pas. C’est 
un écoulement laminaire. 

v Si la vitesse est constante, le mouvement est rectiligne uniforme. Dans ce cas, F est la 


force de frottement. 





Figure I-3 :La répartition des vitesses de fluide placée entre deux plaques 


Le coefficient de proportionnalité u est appelé viscosité dynamique du fluide et V/e 
représente le gradient de vitesse. En effet, une particule de fluide en contact avec une 
paroi prend la vitesse de la paroi. 
Vitesse d'une particule de fluide en contact avec la paroi mobile : V=Ve (y=e) 
Vitesse de la particule en contact avec la paroi immobile : V,=0 (y=0) 

Comme la répartition des vitesses est linéaire(Figure 1-3), la vitesse est proportionnelle à y : 


v=a.y 





eta = — = = 


dV Ve-Vo U 
dy e—0 e 


; s dv 
Cette grandeur est le gradient de vitesse : — 


dy 

2 SV 

Nous avons vu précédemment que :F = u = 
Sdv 
>F = H 
dy 


. A \ sir F dv . 
Si on ramène à l’unité de surface : T = <H a Loi de Newton 


Remarque 
v Ret R;= 1 à quelques centimètres et entrefer R:-R, = quelques millimètres. 


cisaillement constant si R;-R<< R,=(R;-R;)/R;<0,15 entrefer étroit ; 
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v Forme conique de la base du cylindre intérieure pour éviter les effets parasites. Il 
existe deux types de rhéomètres de Couette: 
v À cisaillement imposé: on impose la vitesse de rotation et on mesure le couple C. 
v À contrainte imposée: on impose le couple et on mesure la vitesse de rotation 
obtenue. 
L'expérience de Couettepermet la mesure du coefficient de viscosité par la mesure du 


couple C : 


Ce 


C=R,Su=avecS = 2TRzh = p = E AT 


(2) 


I-2Dimensions et unités 
a-Viscosité dynamique 
De la relationg = = on en déduit la dimension de la viscosité dynamique. 


MTL. K de 
[u] = ara MLL 


En Système International(SI), l'unité de la viscosité dynamiqueest lePa.s ou poiseuille 
(symbole PI) 1 PI1=1 kg.m .s -On trouve aussi le coefficient de viscosité dans le 
Système d’unités (CGS ), l'unité est le Poise(symbole Po) tel que 1PI-10 Po = 1daPo=1 
kgms" 
Autres unités : La viscosité de produits industriels (huile en particulier) est exprimée 
au moyen d’unités empiriques : degré ENGLER en Europe, degré Redwood en Angleterre, 
degré Saybolt aux USA. 

Par rapport aux faits expérimentaux, on est conduit à considérer deux types de fluides : 

v D’une part les fluides newtoniens qui satisfont à la loi de Newton. Ces fluides ont 
un coefficient de viscosité indépendant du gradient de vitesse. C’est le cas des 
liquides purs (eau, huiles). 

v D'autre part les fluides non-newtoniens. Ce sont les solutions de polymères, les 
purées, les gels, les boues, le sang, la plupart des peintures, etc... L’étude de ces 


fluides relève de la rhéologie. 
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b-Viscosité cinématique 


Dans de nombreuses formules apparait le rapport de la viscosité dynamique n et de la 


masse volumique p. Ce rapport est appelé viscosité cinématique v. 


Dimensions : [v]= LT 


=, 
v= 0) 


f , E. l zy 2 -1 
Unités : Système international (ST), l’unité de la viscosité est m .s . 


En Système CGS, l’unité est le Stockes( symbole : St) ms =10" st 


Quelques valeurs typiques 


Quelques valeurs caractéristiques de viscosités dynamiques et cinématiques sont récapitulées 


dans le tableau I-1. 


Tableau 1-1 : Valeurs caractéristiques de viscosités dynamiques et cinématiques obtenues à 


20°C et 1 bar. 


Viscosité dynamiques (Pa.s) | Viscosité cinématiquesv (m’/s) 


H2 
Air 
Hexan 
Benzène 
Eau 
Mercure 
Lait 
Sang de l’homme à 37°C 
Huile d’olive 
Glycérol 
Gels /crèmes 


Vernis /peintures 


Résines/goudrons/bitumes 


Glaces à 0°C 


Granit 
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8,9 x 10° 105 x 10° 
18,5x 10% 15,6x 10% 
0,3x 10° 0,46x 10° 

0,652x 10° 0,741 x 10% 
1,005 x 10° 1,007 x 10% 
1,554 x 10° 0,1147 x 10% 
2,0 x 10° 1,93 x 10° 
4,0 x 10” 
84 x 10° 91,5 x 10° 
1,49 1182 x 10° 
1 à 100 
10 à 1000 
100 à 100 000 
10! 
10% 
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I-3Variation de la viscosité dynamique avec la température et la pression 
Cas des gaz 
Influence de la pression = 0, 499pcl(4) 
oùlest le libre parcours moyen, 
pla masse volumique et c la vitesse moyenne des molécules. 
En général pour les gaz, la pression influe peu sur la viscosité. 
Influence de la température 


Soit uọ la viscosité à la température To 





u _ KNT 

T C 

Ib = 
HT 


oùK et C sont des constantes qui dépendent de la nature du gaz 
Cas des liquides 


Influence de la pression 


pr [KH a 
— = AVEC up elua, les viscosités aux pressions P et Ph 
Pour les huiles minérales. a= 1.003 
Influence de la température 
Toa 
H =r T logii avec ji la viscosité à la température Fy; 


I-4Lois de comportement 


L'expérience montre que, lors de l’écoulement d’un fluide, la pressionne suffit pas à 
expliquer les phénomènes et qu’il convientd’introduire des forces tangentielles qui 
s'opposent au mouvementdu fluide(Figure I-4). Ces forces, de type frottement, dues 
auxinteractions entre molécules du fluide, sont appelées forces deviscosité. La contrainte 
(force par unité de surface) oqu’exerceune couche de fluide supérieure sur un élément de 


surface d’une couche de fluide inférieure,s’écrit : 


—_ dF —_ => 
O = g Inn + 0,taveco, = —p 
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Figure I-4:Les contraintes tangentielles et normales entre deux couches successives de fluide 


en écoulement unidimensionnel. 


a-Fluide newtonien 

Entre deux couches successives de fluide en écoulement unidimensionnel à la vitesse v, il 
existe des contraintes tangentielles à l’écoulement qui accélèrent la couche la plus lente et 
ralentissent la couche la plus rapide. Par définition d’un fluide newtonien, les forces 


visqueuses sont proportionnelles à la différence de vitesse c’est-à-dire au gradient de vitesse: 


Ov s 
o, = pz; Hr) 


„3v E ; ; : z E i 
où désigne le gradient de vitesse dans la direction normale à la surface. De manière 


générale, la contrainte visqueuse varie comme la vitesse de cisaillement; y. 
u est la constante de proportionnalité. Il est caractéristique du fluide et désigne la 


viscosité dynamique du fluide. 


b- Fluides non newtoniens 


Pour les fluides non newtoniens, on peut définir de la même façon une relation de la forme: 
_ pin g 
0, = KŒ = Ky(7) 


OÙ Happ: représente une viscosité apparente. 


La viscosité apparente du fluide n’est pas constante contrairement au cas d’un 


fluidenewtonien, elle dépend du taux de cisaillement. 
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I-5 Relations entre contrainte et taux de cisaillement pour différents fluides 


(Loi de puissance ou d’Ostwald — de Waehle) 


Équation constitutive du fluide dite de type loi de puissance ou d’Ostwald — de Waehlle: 
Happ = ký” "20: = ký”(8) 

oùK la consistance (Pa. s" 

n indice d'écoulement 


v loi valide pour une certaine gamme de taux de cisaillement. ! 


Fluide à seuil de 


contrainte (shear stress) Koarer g 
1 ou plastique 


Gt Fluide pseudoplastique 


(n<1) on rhéofluidifiant, pl 
… (viscoplastic shear-thinning 
Fluide newtonie 


(n=1) (linear newtonian) 


Fluide dilatant 
(n>1) ou rhéoépaississant, T 


(viscoplastic shear-thickening) 





E du/dy taux de cisaillement 
(shear strain rate) 





Figure LS :Les rhéogrammes des fluides. 


La courbeo, = f (appelée "rhéogramme" (Figure 1-5), présente des allures différentes 


suivant selon la nature des fluides. 
a- Fluide newtonien (n=1): 


La viscosité apparente est constante en fonction du taux de cisaillement. C'est le cas de 


nombreux liquides purs, tels que l'eau ou le glycérol. 
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b- Fluides pseudoplastiques ou rhéofluidifiants (0<n<1) 


Les fluides pseudoplastiquessont tels que la viscositéapparente diminue lorsque la 
déformation augmente. Ondit que ces fluides, qui sont les plus fréquents et les plusimportants 
en Génie des Procédés, présentent unerhéofluidification. Appartiennent à cette classe les 
huiles et lesgraisses de lubrification, le pétrole, le sang, certaines peintures, la plupart des 


solutions de polymères. 

c- Fluides dilatants ou rhéoépaississants (n>1) 
La viscosité apparente augmente lorsque le gradient de vitesse augmentetel que les 
suspensions fortement concentrées dans lesquelles la phase liquide n’occupe 
pratiquementque les interstices entre particules solides, les solutions colloïdales d’argile dans 
l’eau, sable mouillé, suspensions concentrées d’amidon. On utilise la loi de puissance pour 
les représenter. 


d- Fluides à seuil ou fluides (visco)plastiques 


Fluide à seuil ou fluides (visco) plastiques tel que lessuspensions,solutions de polymères, 


crèmes et émulsions, pâtes dentifrices, ciment, frais, concentrés de tomate). 


e- Fluides à seuil de type Bingham 


Il existe une contrainte seuil au-dessous de laquelle il n’y a pas d’écoulementtel que le 


dentifrice. 
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II-1 Écoulements laminaires et écoulementsturbulents 
Etudier l’écoulement d’un fluide réel revient à résoudre l'équation fondamentale de la 
dynamique (équation de Navier-Stokes): 

dv = — — 
a TVP + HAD + pg0) 
Mais en pratique, cette équation ne peut se résoudre analytiquement qu’en introduisant des 
hypothèses de travail qui permettront de résoudre cette équation dans le cadre de régimes 
d'écoulement particuliers :le régime laminaireet le régime turbulent. 
v On dit qu’un écoulement est laminairelorsque le mouvement des particules fluides se 
fait de façonrégulière et ordonnée. 
v  L’écoulement est turbulentlorsque le déplacement est irrégulier et que 

desfluctuations aléatoires de vitesse se superposent au mouvement moyen du fluide. 

En régime turbulent, on devra utiliser des relations empiriques généralement 

déterminées expérimentalement. 
L'une des premières analyses de la transition d’un régime laminaire vers la turbulence est 
basée surdes observations d’écoulements dans une conduite cylindrique effectuées par O. 
Reynolds en 1883. 
Le montage expérimental schématisé sur la figure H-2 comporte un réservoir de liquide sous 
pression débouchant sur un conduitcylindrique. Un tube mince permet l’injection de colorant. 
Lorsque l’écoulement est laminaire, le filet decolorant reste mince, régulier et parallèle à la 
paroi du cylindre. En écoulement turbulent, le colorant estrapidement dispersé. Dans cette 
situation, une mesure de la composante de vitesse axiale montre quecelle-ci fluctue de façon 


aléatoire dans l’espace et dans le temps. 





Figure II-1 : Expérienced'OsborneReynolds 
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laminaire 6 + turbulent 


(v (t) et v (t) << v_{t) 


en un point M : < 
v 


Vz 


laminaire 
turbulent 





Figure I-2 : Différents régimes d’un écoulement 


À partir de ces observations, les travaux de Reynolds ont alors permis de montrer que la 
transition du régime laminaire au régime turbulent (Figure II-2)n'est pas seulement 
conditionnée par le débit mais dépend d'un ensemble de paramètres qu'il convient de 
regrouper pour former une quantité sans dimension qu'on appellera «nombre de 
Reynolds »Eq. (10). Ce nombre prend en compte la vitesse moyenne de l’écoulementV, le 
diamètre Dde la conduite (ou tout autre paramètre de longueur caractérisant l'écoulement), 


ainsi que les propriétés intrinsèques du fluide (masse volumique pet viscosité p) : 


Re == (10) 


où n: masse volumique [p ] = ML” 
u : viscosité [u ] = ML” T 

v= w/p : viscosité cinématique [ v] = L °T” 
V : vitesse [v ] = LT” 
D : diamètre [D] = L 
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Remarque :Compte tenu des dimensions de chacun des paramètres, on vérifiera facilement 


que le nombre de Reynolds est bien sans dimension. 


On constate généralement que la transition d’un régime laminaire à un régime 
turbulent(Figure I1-3) 
s’effectue lorsque Re = 2000 = Re., nombre de Reynolds critique. 
v Pour Re< 2000, l’écoulement reste laminaire et une perturbation localisée introduite 
dans l’écoulement est progressivement dissipée. 
v Dans un intervalle de Re de 2000 à 3000, des « paquets » turbulents sont convectés 
dans le conduit de façon intermittente. 
v Aux nombres de Reynolds plus élevésRe> 3000, l’écoulement devient turbulent dans 
son ensemble ; c’est-à-dire que les forces de viscosité ne sont plus suffisantes pour 
empêcher les inévitables perturbations d’engendrer des tourbillons qui se superposent 


à l’écoulement global. 





Figure I1-3 :Les régimes d’écoulement 


Remarque : 

En prenant certaines précautions pour éviter au maximum les perturbations (vibrations, 
aspérités de la conduite...) l'écoulement peut rester laminaire au-delà deRe,, jusqu'à une 
certaine limite (imprécise) qui dépendra essentiellement du niveau de précautions prises pour 


repousser l'apparition des turbulences. 


II-2 Écoulement laminaire et pertes de charge régulières 
Que l'écoulement soit laminaire ou turbulent, la résolution de l'équation de Navier-Stokes 


passe par l'introduction de la notion de pertes de charge. Ainsi, pour rendre compte de la 


| 11 | 
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dissipation d'énergie due aux frottements visqueux, ces pertes de charges prendront place 
dans la formulation d'une équation deBernoulligénéralisée. 
Partons de l'équation de Navier-Stokes(Eq. (9)) établie pour l'écoulement d'un fluide 


newtonien incompressible : 
= az Ç . OÙ T | : S S avr 
Si l'écoulement est stationnaire = 0. la dérivée particulaire se résume à la dérivée 


convective et peut se développer comme: : 


0 


EA 


dv œ vo = lv Y 
S -= + @V)5 = (VV )P = < VGS -VA VAD) Fe 


dt 
On peut alors écrireEq. (12) 
1 — Dj rm, 2 = G 
zPY — pÜA(VAT)= —Vp + uA? + pg (12) 
L'accélération de la pesanteur dérivant d'un potentiel et le fluide étant incompressible (p= 
Cte), on regroupe sous un même gradient trois des cinq termes de cette équation pour lui 


donner la forme suivante Eqs. (13-15) : 
—/1 EER E 2 cs ni. 
V (or) — pDA(VAT)= —Vp + uA? + VpgZ (13) 
e: a 1 = > r = 
V(p+ 092 +3pv?) = HAD L 00 AUV AY] (14) 
ca X 1 m E 
V (p + pgZ + zpr?) = pA? + 2pv A Q(15) 
On retrouve ainsi une formulation semblable à celle conduisant pour un fluide parfait à 
l'équation de Bernoulli. 
De la même manière, le produit vectoriel du vecteur vitesse avec le rotationnel de la vitesse 


s'annule lorsque l'équation est considérée sur l'ensemble des points constituant une ligne de 


courantEq. (16). 


EN n 1 = = lan 6 
V (p + pgZ + zpr?) dse, = (HAV + 2pv ^ Q) dse, (16) 





(510) 1 Tete, || Ð = (74019, = 0 
Il reste donc que le long d'une ligne de courant, l'écoulement stationnaire d'un fluide de 


viscosité non négligeable obéit à l'équation suivanteEq.(17) : 


= E 1 2 S 
V (p + pgZ + zpr?) dse, = uAvdse, (17) 
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Puis la projection sur chacun des axes d’un repère cartésienEqs. (18-23) : 
dse, = dye, +L dye, + dze,(18) 


2C)+ dtd 
sB + gy &(19) 











(e) B 
V(e) = zy x + 





VU = VD-6- + VV e, + VV,e,(20) 


d 1 
S a = pv? — 
3x (p +pgZ + 2 pv ) dx = uAv,dx(21) 


ð — 1 
F (p + pgZ + pv?) dy = uAv,dy(22) 
Zp + p9Z +Żpv?)dz = uAv,dz(23) 


SoitEqs. (24-26) après simplification pardx, dy et dz : 


ð EE: 
2 (P +p9Z +5 pv ) = HAD, (24) 


d z 1 
a (p + pgZ + pv?) = uAv, (25) 
d = 1 
2 (p+pgZ +=pv?) = uAv,06 
p + pgZ + = pv?est la pression totale (également appelée « charge »), notée P;,, somme 


de la pression motrice p + pgZ et de la pression cinétique = pv? . C'est donc une quantité 
scalaire qui doit vérifier l'équationEq. (27): 

Vp, = HAT (27) 
On suppose alors qu'un repère cartésien soit choisi de telle sorte qu'en tout point de 
l'écoulement on ait v=u(x,y,z)e, et que l’écoulement est laminaire. Dans ces 
conditions, la projection dans les trois directions du repère de l'équationEq. (27) 


précédemment établie, donne : 


Vv, =v 
y=0>V=vVve, 
v,=0 
ôPt 
Jı T KAV 
ôp: 
dy 
dH: 
—— = 0 
0z 
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On en déduit que la pression totale n'est fonction que de x, et qu'elle varie en obéissant à 


l'équationEq. (28) : 
‘dx _ nL T ðy? t 71 (8) 


Par ailleurs, l'écoulement étant conservatif, l'équation de continuité Eq. (29) s'applique dans 


sa forme la plus simple: 


== ðu ðv ðw ðu 
Vv=0>— + —+ —=0 > — =029 
ôx + dy T oz ox ( ) 
: r ðu + 12 
et commev=w=0, il en découle Pi 0 >u=u (y,z). Nous pouvons alors reconsidérer 


l'équation régissant l'évolution de la pression totale en la formulantEq. (30) : 


dp: _ 


3u u 
FF = u + = CteV(x, y,z)(30) 


dy? 077 
Il s'avère alors que le membre de gauche E5 ne peut dépendre que de xalors que celui de 


. 92 92 
droite (+ = 


S 2 Jne peut dépendre que de y et z. La conséquence est donc que les deux 
membres sont égaux à une même constanteW(x, y, Z). Ceci permet de montrer que la pression 
totale Eq. (31) varie linéairement avec x: 


dp, 
— = Ct 0 31 
dx e< (31) 


Remarque : 
v Pour un fluide parfait dont la viscosité est négligée, la pression totale ne varie pas. En 
effet, la constante étant proportionnelle à u, elle devient nulle, l'équation de Bernoulli 


Eq. (32) s'applique le long d'une ligne de courant : 


SL Z 0>p+pgZ +5pv°= Cte(32) 


On dit alors que la charge est constante ou bien encore qu'il n'y a pas de perte de charge. 

v Pour un fluide réel (dont la viscosité n'est pas négligeable), l'écoulement laminaire est 
caractérisé par une variation linéaire de la charge avec la distance parcourue. En effet, 
puisque les frottements visqueux sont responsables d'une dissipation d'énergie et que la 
pression totale n'est autre qu'une énergie volumique, il s'ensuit logiquement que la 


pression totale doit diminuer avec la progression de l'écoulement. 
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Figure II-4 : Evolution de la pression totale (charge) avec la progression de l'écoulement 


Dans l'hypothèse logique d'une pression totale diminuant linéairement avec la progression de 
l'écoulement, le graphe de la figure H-4 montre que sur une longueurAx=x;-x, la chute de 


pression totaleAP est proportionnelle à Axet s'exprimeEq. (33) : 


dp, Pe- Pa _ AP: 
dx X2 — X: Ax 


> Ap, = — ŽŽ Ax63) 


Ainsi, une chute de perte de charge se traduit parAp, > 0 et est la conséquence d'un gradient 
de pression constant et négatif <0. C'est pourquoi on parle de pertes de charges 
régulièresEq. (34): 


d 
> Pu = Pe- mE (x2 — x1) = Pe TAH > U (34) 


On peut généraliser l'équation deBernoulli en y faisant apparaître les pertes de charges 


régulières de la manière suivante Eq. (35) : 


= 1 — 1 
Pı + P9Z1 +5 PV1° = P2 +P9Z2 +3PV:2* + Api(35) 


oùAp, représente l’ensemble des pertes de charge entre (1) et (2). 
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II-3 Coefficient de perte de charge 
Particulièrement, les pertes de charge entre l'entrée et la sortie d’une canalisation est la 
somme de deux types de contributions : 
v les pertes de charge systématiques (appelées aussi linéaires ourégulières) dues aux 
frottements visqueux. 
v les pertes de chargesingularités dues aux changements de direction et de sections 


(rétrécissements et élargissements) dans la conduite; 


I1-3-1 Perte de charges linéaires ou régulières 
Les pertes de chargesystématiques (linéaires ou régulières) (chute de pression Ap = pı - p2) 
résultent du frottement exercé entre le fluide et la surface intérieure de la canalisation. Elles 
sont proportionnelles à la longueur L de la conduite et au carré de la vitesse moyenne V du 
fluide, inversement proportionnelle au diamètre d et fonction de la rugosité moyenne e de la 
canalisation. Entre deux points séparés par une longueur L, dans un tuyau de diamètre D 
apparaît une perte de pression Ap. Exprimée sous la forme suivante Eq. (36) : 
Ap = 170 Lae 

où D: diamètre de la conduite considérée (m) 

L : longueur de la conduite considérée (m) 

V : vitesse moyenne (m/s) 

(Sans dimension) : coefficient de perte de charge linéaire. Il dépend de la nature de 

l’écoulement et de l’état de surface de la conduite. 

Remarque : 

La nature de l’écoulement est caractérisée par le nombre de Reynolds Re 


v l’état de surface est défini par l’épaisseur moyenne des rugosités : 
v Alors f (Re, g/d) 


Expérience de Nikuradze 


L’expérience de Nikuradze(Figure Il-5) a pour but de déterminer l’influence de la rugosité 


des parois sur le coefficient de perte de charge linéaire. Les parois d’une conduite sont 


| 16 | 
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rendues artificiellement rugueuses en y collant des grains de sable calibrés. En changeant la 
taille des grains on change la rugosité et on trace log à enfonction de log Re. 
On appelle ks la taille moyenne des grains collés et D le diamètre de la conduite. g /D est la 
rugosité relative. On remarque sur les abaques de Nikuradze(Figure I-6 ) que: 
v la rugosité n’a pas d’influence en régime laminaire et le nombre de Reynolds critique 
RC est indépendant de la rugosité. 
v quand Re>Recla perte de charge dépend fortement de e /D ou du nombre de Reynolds 
v quand R >>>Re la perte de charge ne dépend pas du nombre de Reynolds, elle ne 
dépend que de £ /D. 


Expérience de Nikuradze 
Reynolds + modification de l'état 
de la canalisation en collant des 
grains de sable. 


réservoir 





pogano EDO G-DBLD opoo 


D (m= H G 





Resistance coefficient, f 











Figure I1-6 : Abaques de Nikuradze 


Selon le nombre de Reynolds, on distingue différents cas : 
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a- Cas de l'écoulement laminaire (Re< 2000 -la loi de Poiseuille) 
Dans ce cas on peut montrer que le coefficient À est uniquement fonction du nombre 
deReynolds Re; l'état de la surface n'intervient pas et donc À ne dépend pas de de la rugosité 


g(noté aussi k), ni de la nature de la tuyauterieEq. (37) . 
64 
À = FASI) 


b- Cas de l'écoulement turbulent (Re> 3000) 
Les phénomènes d'écoulement sont beaucoup plus complexes et la détermination 
ducoefficient de perte de charge résulte de mesures expérimentales. C'est ce qui explique 
ladiversité des formules anciennes qui ont été proposées pour sa détermination.En régime 
turbulent l'état de la surface devient sensible et son influence est d'autant plus grandeque le 
nombre de Reynolds Re est grand. Tous les travaux ont montré l'influence de la rugosité eton 
s'est attaché par la suite à chercher la variation du coefficient à en fonction du nombre de 


Reynolds Reet de la rugosité. 


La formule de ColebrookEq. (38) est actuellement considérée comme celle qui traduit le 
mieux lesphénomènes d'écoulement en régime turbulent. Elle est présentée sous la forme 


suivante : 


1 E 2,51 
Ti = -2log Gp + Reya 68) 





L'utilisation directe de cette formule demanderait, du fait de sa forme implicite, un calcul par 
approximations successives ; on emploie aussi en pratique des représentations 
graphiques(abaques). 

Pour simplifier la relation précédente, on peut chercher à savoir si l'écoulement 
esthydrauliquement lisse ou rugueux pour évaluer la prédominance des deux termes 


entreparenthèses dans la relation de Colebrook. 
Remarque : 


On fait souvent appel à des formules empiriques plus simples valables pour des cas 


particuliers et dans un certain domaine du nombre de Reynolds, par exemple : 
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a- Formule de Blasius Eq. (39) : pour l’écoulement turbulenthydrauliquementlisse et 
3000<Re< 10”,Àne dépend que de Re : 
À = (100Re) °?5 = 0,316.Re °?°(39) 


b- Formule de Blench Eq. (40): pour l’écoulement est dit turbulent 


hydrauliquementrugueux et Re>105, Ane dépend que de s/d . 


À = 0,79 [Eco 


oùs : rugosité de la surface interne de la conduite (mm) ; 

d : diamètre intérieur de la conduite (mm). 

Diagramme de Moody 

Cette relation est appelée équation deDarcy-Weïisbach. Le facteur adimensionnéAdésigne le 
coefficient de perte de charge régulière et ne dépend que du nombre de Reynolds pour une 
canalisation lisse. Dans le cas d’une canalisation rugueuse, un quatrième nombre sans 
dimension intervient : la rugosité relative #/Dqui mesure le rapport de la hauteur moyenne 
des aspérités de la paroi interne de la conduite sur son diamètre interne. La valeur de Apeut 


être obtenue à l’aide d’abaque comme le diagramme de Moody(Figure I1-7). 





0.05 
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Figure II-7 : Diagramme de Moody 
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IL-3-2 Perte de charge singulière 

Les pertes de charges singulières résultent de la présence de coudes, raccords, 
branchements, robinets, etc. Tous ces éléments (singularités), installés le long des 
canalisations, constituent des obstacles qui freinent le passage du fluide et amènent des pertes 
de charge. Les pertes de charge singulières sont proportionnelles au carré de la vitesse, elle 


est exprimée sous la forme suivante Eq. (41) : 


y2 
Ap = Kp 749) 


où K :Le coefficient de perte de charge singulière (sans dimension). Il dépend de la 
nature de lasingularité. 
v : La vitesse moyenne incidente du fluide arrivant sur l’obstacle. 
Remarque : 
La plupart des singularités ne permettent pas une approche analytique simple comme celle 
qui vient d'être développée concernant l'élargissement brusque. Les coefficients de perte de 
charge sont donc le plus souvent déduit de mesures ou de calculs numériques, et sont 
disponibles sous forme d'abaques(Figure II-8). À titre d'exemple, la figure H-9 liste de 
façon non exhaustive quelques singularités typiques avec les coefficients de perte de charge 


associés. 


Conduite hyciauliqguement 


hydrauliquement lisse 





Figure IT-8 :Abaque des pertes de charge. 
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E 
Elargissement brusque Divergent 
K =K1-S, S, X =-Q0-S, SŸ sin 





Rétrècissement brusque Convargent 
K <U -i}° K -Ü 01T sin 
F H S L P 1 


ug 


Coude brusque Coude arrondi 


K = sin ter 1 Zain" © 


A N [a31 . L842LD RY >] 








Entrèe de canalisation brusque Entrée de canalisation progressive 
K - 05 K = 0,0 


[= = 


Figure IT-9 :Liste non exhaustive de quelques singularités typiques 





IL-3-3 Pertes de charge totales 


Pour un réseau hydraulique complet comportant différents tronçons de longueurs et sections 

différentes et reliés entre eux par des singularités, on peut évaluer les pertes de charge totales 

entre l'entrée et la sortie du circuit en formulant l'équationdeBernoulligénéralisée Eq. (42) : 
Li 


> 4 ; 1 
AP: = SHH Air + Dit > PV” Ki(42) 
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HI-1 Introduction 

Par rapport aux écoulements externes, les écoulements internes ont ceci de particulier que les 
conditions aux limites y sont omniprésentes. Ils sont abordés dans ce chapitre sous deux 
aspects complémentaires : les profils de vitesse et les pertes de charge qui conditionnent le 
dimensionnement des circuits fluides. 

Les écoulements internes ont pour caractère particulier l’absence de conditions aux limites à 
linfini, ce qui n’est pas le cas des écoulements externes. L’influence des parois se manifeste 
donc dans toutes les directions, sauf éventuellement aux sections d’entrée et de sortie du 
fluide, et en particulier le dessin deséchangeurs ; plus généralement, les dissipations doivent 


être prises en compte dans les bilansénergétiques ou économiques. 


II-2Ecoulement parallèle dans une conduitebidimensionnelle 
On considère une conduite bidimensionnelle infinie, d’épaisseur e=2b faible, l’axe x coïncide 
avec l’axe de symétrie de la conduite. Le fluide est considéré comme newtonien, 


incompressible en écoulement stationnaire (permanent) selon l’axe x(Figure III-1). 





Figure III-1 :Écoulement parallèle dans une conduite bidimensionnelle 


On choisit un repère cartésien, la projection de l’équation de Navier Stokesaprès 


modification, donne Eqs. (43-45) : 
ð 1 
ka TZ +- pv? |= 
x (p +pgZ +5PY ) HAD, (43) 
b) =à 1 
(+ 007 +4pv?)= uav, 


d = 1 
T (p +pgZ + pv?) = uAv,(45) 
Dans ces conditions, la projection dans les trois directions du repère cartésien de l'équation 


précédemment établie, donne : 


[a] 
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Vv, D 


D, = 0 
leséquation précédentes deviennent : 
dpt _ 
Fe = uAv 
dpi _ . 
0y 
ôP: 
>t- 0 
oz 


On en déduit que la pression totale n'est fonction que de x, et qu'elle varie en obéissant à 
l'équation Eq. (46) : 


dp: 3u ðu ðu 
dx ðx? dy? 8z? 


(46) 
v Symétrie : dans tout plan xOy, on a le même profil donc w = 0 
v Ilny apas de composantes de la vitesse suivant y donc v = 0 


v L’équation de continuité div V = 0 impose du/dx=0 


Le profil des vitesses dans la conduite est donc donné par la résolution de l’équation 
différentielleEq. (47) : 

“Pr = ga = Cte (47) 
On intègre deux fois cette équation et on obtient Eq. (48): 
2 
u 


u(y) = + (2) E Ky + K'(48) 


dx 

oùK et K’ sont des constantes d'intégration, elles sont déterminées grâce aux conditions 
aux limites sur la vitesse ( une particule de fluide en contact avec une paroi possède la 
vitesse de la paroi) soit:u(b) = u(-b) = 0 

En insérant ces deux conditions dans l équation précédente, on obtient un système de deux 


équations à deux inconnues facile à résoudre, on obtient Eq. (49): 


= r _ 1 (dPt\p2 
K=0 et K = IGL (49) 
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En reportant dans l’équationEq. (48), on trouve le profil des vitesses suivantEq. (50) : 
1 rd 
uy) = ZL] 0 - Eo 


Le profil des vitesses est parabolique. 





Figure II-2 :Le profil parabolique des vitesses de fluide en écoulement parallèle dans une 


conduite bidimensionnelle 


Le gradient de pression peut être évalué, en effet si on néglige les forces de pesanteur Eqs. 


(51-52): 
Da P = P(x) 51 
ay de. P = P SD 
_— dp 


On peut donc écrire en appelant P, et P les pressions à l’entrée et à la sortie de la conduite 
Eq. (53): 


dP Pı — P2 _P17 P2 








(53) 


dx Xu Xa L 
L est la longueur de la conduite et P,>P,; 
Il suffit de mesurer les pressions correspondantes pour déterminer le gradient de pression 


La vitesse moyenne de l'écoulement est définie parEq. (54) : 


n< LH US 
H< U VdS (54) 


oùS : La section de la conduite(m”)etqy : débit en volume (m°/s) 
Le débit en masse par unité d’aire Eq. (55): 


Ka = pU(55) 


Pour le conduit bidimensionnel, on trouve : 


La vitesse maximaleEq. (56): 
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— _ 1 fdp:\ p2 
Umar = -= (2) b (56) 


La vitesse moyenneEq. (57): 


2 


e (2) b? = Šumas (57) 


Le débit en masse par unité d’aireEq. (58): 


S = g OT Ob (68) 
III-3Écoulement de Couette 

Une autre situation simple est celle de l’écoulement de Couette(Figure IMI-3). On a les 
mêmeshypothèses que précédemment sur la géométrie de la conduite, sur le fluide 
etl’écoulement. La différence vient du fait que la paroi supérieure possède unevitesse Vs. 
Comme dans l’exemple précédent, on néglige les forces de pesanteur. 

On obtient bien sûr la même équation différentielle avec des conditions auxlimites qui 


changent. 





Figure III-3 : Écoulement de Couette 


Les conditions aux limites deviennent : V(+b) = Ya et V(-b) = 0 


La solution générale de l’équation est comme dans l’exemple précédentEq. (59) : 


2 
V(y) = = (2) Z + Ky + K'69) 


Avec les conditions aux limites, on trouve K et K’Eq. (60) 


= Vo r= Vo 1 (dr) p2 
— 2b A ($b 199) 
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La solution s’écrit doncEq. (61): 


1 


u(y) < URT -y + 2 Y+b) 6D 


| | 


(Poiseuille)(Couette) 


On peut écrire le profil des vitesses sous la forme : 


On pose : 


k= (8) v < B+ D + ne 


Le profil de vitesse est donc parabolique. 
v Kiapparaît comme un gradient de pression adimensionnel. 
v K;> 0 le gradient de pression décroissant est dans le sens du mouvement doncu(y) est 
toujours positif. 
v K< 0 et inférieur à —1/2 la vitesse devient négative au voisinage de la paroi(le sens de 
l’écoulement est inversé). 
En l’absence d’un gradient de pression (dp4/dx = 0) la répartition de vitesse entre les deux 


plans devient linéaire Eq. (63): 
V(y) = 7 (y +b) sachant que e = 2b(63) 
Le cisaillement ©; dans l’écoulement est alors uniformément répartiEq. (64) : 
O, = n = p” = Cte(64) 


Le débit total de fluideentre les deux plans est dû à la fois à la vitesse d’entraînement et au 
gradient de pression. En désignant par /la largeur de la veine fluide, le débit-volume q, 


s’écrit, dans une section SEq. (65): 


qy = f uds = l ff udy(65) 
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On obtient en intégrant et en groupant les termes Eq. (66): 


e e? 
qv = l5 Ue — ap Pert) (66) 


III-4Écoulement de Poiseuille 


On considère un écoulement horizontal, laminaire, incompressible et newtonien d’un fluide 
de viscosité pm, de masse volumique p, dans un tube cylindrique de rayon R 


petit.L’ écoulement est permanent dans la direction x(Figure IIT-4). 





Figure III-4 :Écoulement de Poiseuille(Écoulement parallèle dans une conduite 


cylindrique) 


La vitesse de la couche de fluide en contact avec les parois est nulle. Elle augmente lorsque 
la distance r à l’axe diminue. 
La vecteur devitesse en tout point à l’intérieur de la conduite s’écrit Eq. (67): 
V(r, 0, x)é,carv,=v.=0(67) 
L'écoulement étant conservatif et le fluide incompressible, l'équation de continuité se résume 


à VY = 0 se développe en coordonnées cylindriques comme Eq. (68): 
13 1 va ðu 
Tir (rv,) + ET + Jx = 0(68) 


Et puisque v,=vọ=0, il reste Eq. (69): 


ðu 
e 0(69) 


traduisant le fait que la vitesse reste constante le long de l'axe de la conduite. 
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us ; à & ee ; ð 
Compte tenu de la symétrie de révolution du problème, on en déduit par ailleurs que = 0, 


ce qui permet d'en conclure que le vecteur vitesse est purement axial et ne dépend que de 
rEq. (70): 

V(r)é,(70) 
D'autre part, l'équation régissant l'évolution de la pression totale (charge) P, le long d'un 


écoulement laminaire orienté suivant l'axex (voir section précédente) Eq. (71): 


dP; _ 
= = pAu(71) 
En cordonnées cylindriquesEq. (72): 

roue ne 
dr r2 002 0x2 rôr 


alorsEq. (73): 


dE _ 10 (p? 
‘dx ré T or (r loa) 


où le terme de gauche ne dépend que dexet celui de droite que de r. 


Il s'en suit que les membres de cette équation sont égaux à une même constante K. Ainsi, on 


a Z = K traduisant la perte de charge régulière ; 
1 0 / ðu C | | 
u RER (r =) = K > permettant la détermination du profil de vitesse u (r). 


Résolvons cette dernière équation pour établir le profil de vitesse sur une section de la 


conduite Eq. (74): 


13 ou 1 0 ðu Kr ðu Kr? ðu Kr K' 
P(r) er LR RA LN aD ELT 
dr H dr 2u or 2u r 


et donc : 


Kr? 
u(r) = An + k'In(r) + k" (74) 


oùk'et k" sont deux constantes d'intégration dont l'indétermination peut être levée en 
utilisant les conditions aux limites. 
Ainsi, on sait que la vitesse des particules fluides est nécessairement nulle au contact avec les 


parois de la conduite. On a doncEq. (75) : 
r2 
u(R =r)=0> T + k'In(r) + k" = 0(75) 


Par ailleurs, la vitesse à l'axe de la conduite, enr=o, soit finie, il faut que K’=0alorsEq. (76): 
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Paen (76) 
et que le profil de vitesse se formule Eq. (77): 
ulr) =- Rr’) où O<r<R (77 ) 
Cette équation traduisant leprofil parabolique de u (r)(Figure IMI-5)et que la vitesse 


maximale sur l'axe estEq. (78) : 
K 


Umax = U (o) = um (78) 


u (r)> Oet K<0 imposeK = 2 < 0 Perte de charge. 


RA 





Figure II-5 :Le profil parabolique des vitesses de fluide en écoulement de Poiseuille 


Pour caractériser les pertes de charges, on définit la vitesse moyenne de l'écoulement 
parEqs. (79-80) : 
qy = f uds = L; udy(79) 


Sao a 1 
T = | = —z [J 0nds = — ff u(r) ds (80) 


avecds=2rr dr d'où Eq. (81): 


R R 
u (r)2rnrd a] R? 2\rd 
= —; | uļlr)2rnrdr = —— TdT 
“= TR l E au J, ) 
— _ 2K RH K 
u = auR? 4 au À (81) 
et donc Eq. (82): 
sr — — K p2 — “max 
QT =e 


Et commegq, = US = urR?, on a par ailleursEq. (83) : 
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K 
qv = — pp TR* (83) 


Cette dernière relation permet ainsi de faire le lien entre le débit volumique et les pertes de 





. H i z dP: z _aPt RERET s qv 
charge régulières puisqueK = q, €t que AP; = PP Ax = —KAxoù K = -8u TRE On 


obtient alors la formule de PoiseuilleEq. (84): 





qv H 
AP, = 843 Ax = Su Ax (84) 
La différence de pression totale (ou perte de charge) engendrée par les frottements visqueux 
dans une conduite de longueur L=Axest donc proportionnelle à L, mais également au débit et 


à la viscosité du fluideEq. (85) : 





AP, = 8u LL (85) 


IHI-SÉcoulement dans un espace annulaire 


Examinons encore le cas d’un écoulement entre deux tuyaux droits coaxiaux, que 
l’onrencontre en particulier dans certains échangeurs, et aussi dans les forages pétroliers. 
Onnotera R;le plus petit rayon et R;le plus grand.La figure IIT-6 illustre un écoulement de 


fluide dans une canalisation annulaire. 





Figure III-6 : Écoulement dans une canalisation annulaire. 
La solution générale reste valable, les conditions aux limites étant ici : 


r=R, ,U=0et r=#R,,U=0 
soitles équationsEqs. (86-88) : 
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AP, = -= L(86) 


_ _ RidPe 
AlnR, +B =- (87) 


ALnR; + B = — PA Pu(s8) 
On en déduit les constantes A et BEqs. (89-90): 
A TE 1 dP: (R3-R) R?) (89) 


4u dx Ln Za 


L dP (R5-RŸ) 2 
B = 14 T pi in G LnR, — rilo 
ce qui donne pour U, après regroupement de certains termesEq. (91) : 


— _ 1a 2 (RŻ- 2 
u(r) = PEE [e -r?° + pE n2 Ln z Jon 
1 


Le maximum de la vitesse est obtenu pour un rayon rmtel queEq. (92) : 


(92) 


Y maw = 





En ce qui concerne le débit dû au gradient dP/dx, les calculs sont un peu pluslaborieux. 
Posons provisoirement, pour alléger l’ écriture Eq. (93): 
(R3 — RŸ) 


C= (93) 
2Ln (PE) 


Le débit-volume s’écrit Eq.(94): 


R2 m dP 
q, = Í U2nrdr = -— St “R — r? + CLnr — CLnR,)rdr (94) 
r1 24 dx Jle 


Au total, on a doncEq.(95) : 


RPR RIR} _ R2 Rî_C(p2_ p2 R2 
b=- I - ae - (R2 Ri)EnR; + C foa rlnrdr}(5) 
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Au total, après avoir remplacé C par (Eq.(93)), puis groupé et réarrangé un peu les termes, on 


a doncEq.(9%6) : 





a dP R3-R? 
ER - RD) fri + ri - iai (96) 
1 
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IV-1 Introduction 


Quant le système étudié est trop complexe pour permettre une résolution complète des 
équations fondamentales, ou bien lorsque son comportement est chaotique, l’analyse 
dimensionnelle donne accès de façon simple à des relations entre les grandeurs caractérisant 
ce système. 

Le regroupement de ces différentes grandeurs en des nombres sans dimensionpermet par 
ailleurs d’établir des similitudes entre les comportements des systèmes semblables mais 
différents (prototype/maquette). 

Prenons l’exemple de la détermination des pertes de charge régulières dans une conduite 


cylindrique. 


= har Aa 


Grandeurs: 


Géométriques : 


Physiques : 





Cinématiques et dynamiques: 


Les différentes grandeurs qui interviennent sont : 


AP a7 
— : Perte de charge par unité de longueur, 


L 
D: Diamètre de la conduite, 

v __:Vitesse moyenne del’écoulement(ou Q: Débit volumique), 
p :Masse volumique du fluide, 

u : Viscosité dynamique du fluide, 

g :Rugosité moyenne de la conduite. 


Par conséquent, il existe une relation entre ces différentes grandeurs Eq. (97): 
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APC 
L 





f(D,v;p;u.E)(07) 

La fonction f peut s’avérer difficile à trouver ; l’analyse dimensionnelle va alors nous 
permettre d’établir une relation plus simple entre un nombre moins important de grandeurs 
sans dimension. 


Une méthode systématique va permettre de trouver trois nombre sans dimensionEq.(98) 


AP D 

17 L pr 
Tz = P = R,(nombre de Reynolds)(98) 
Ti = S = £,(rugosité relative) 


On peut aussi établirEq. (99) : 


T1 = $ (T2,T3) 


AP; D vD z AP 
t l 


>—— = Ÿ u D = $ (Re £+) (99) 


L’analyse dimensionnelle permet aussi de voir que la perte de charge régulière est fonction 
uniquement de nombre de Reynolds et de la rugosité relative de la conduite. 
Voyons maintenant une description détaillée de la méthode. 


IV-2 Théorème de Buckingham Théorème de TT ) 


Si une équation comportant k variables est homogène, elle peut être réduite à une relation 
entre (k-r) produit indépendant sans dimension, ou r est le nombre minimal de dimension 
requis pour décrire les k variables. 


Inventaire des variables 
L’inventaire des variables est regroupé dans le tableau IV-1. 


Tableau IV-1 : Inventaire des variables 


Symbole Dimension (MLT) 
Perte de charge 


Diamètre 
Longueur 





L 


Elaziouti Abdelkader 


UMT/FT/ GP/1*° Année Master: Gestion des Déchets Module :Dynamique des fluides réels 


L Viscosité 1" m "| ML T 
L Mase volumique H ML” 


Afin d'illustrer cet énoncée, reprenant l’exemple d’introduction. 


E l AP ns l ai l 
Nous avions k=6 variables Cu D,v,p,u,e),qui nécessitent un minimum de r=3 (trois 


dimensions (M, L, T). 


Par conséquence, l’équation reliant les 6 variables peut êtreramener à une équation reliant 


k-r=6-3=3 produits sans dimension Eq.(100): 


T, = C = R,(100) 


Es. 
= Er 


Le théorème de Buckingham x, permet donc le passage Eq.(101):: 


APC 
L 





f(D,v,p,1,E)=7T 1 = d (T2,T3) (101) 


Pour effectuer une analyse dimensionnelle, on doit considérer les étapes suivantes: 


6- 


Dresser la liste de toutes les variables du problème=k, 

Ecrire l’équation aux dimensions de chacun des K variable, 

Déterminerr, et donc k-r—= le nombre de produit sans dimension caractérisant le 
problème, 

Parmi les k variables,on choisir un nombre r, qui soient dimensionnellement 
indépendant=r variables primaires, 

Former les k-r produits (x) en combinaison les k-rvariable non primaires avec les r 
primaires. 


Formuler la relation entre les k-r produits (x) trouvé. 


A titre d’illustration, appliquant la méthode à un exemple concret.On appelle force de 


trainée, la force Dexercée par un écoulement sur un objet, dans la direction parallèle à 


l’écoulement. Nous allons étudier le cas d’une plaque plane rectangulaire. 
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1. Les variables du problème sont : D, h, L, v, p, =k=6 ; 
D :Force de trainée, 
h : Hauteur de la plaque, 
L : Largeur de la plaque, 
v : Vitesse moyenne de l’écoulement, 
p : Viscosité de fluide, 


u :Masse volumique de fluide. 


2. Equations aux dimensions 
[D] :MLT° 


[h] :L 

IL] :L 

[v] : LT! 

[pl :ML° 

[u] :ML'T' 
3. Nombre de produit x sans dimension 


k-r=6-3=3 


4. Chois de r=3 variable primaire dimensionnellement indépendantes. Par exemple h, v, p. 


Remarque :On ne peut pas choisir à la fois h et L 


5. Formuler les 3 produits (7). 


Par combinaison des variables primaires et non primairesEq.(102): 
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Tı = Dn oP ori 
Tt = Dh% pf? vY? (102) 


Ta = Dn p 393 


Tı = Dh“ oB ort. sans dimension 
0=1+ß; Bi=-1 





D 
=] + .— + = = —— 
> 0=1 (LA 3PB: yı > (oSI 2 TT] BT 


0=—2—y 1Y 1=-2 


ra = Dh®? pf?vY?— sans dimension 
0=ß2ß2=0 


Bá 0=1+0—-3f2+y > o=-1 DT) = z 


0=—2y 1 Y2=0 


ra = Dh%pl3v8-, sans dimensions; = SR 


6. Formuler la relation entre les 3 produits (n)trouvés : 


AP 
G f(D, V, P, H, EIT, = p (T2, T3) 





AVECTT1 = 2 
L 
H = h 
re 
3 pvh 


soit Eq.(103)D = pv'h”Œ(L/h,1/ Re) 
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(Facteur de forme)(Nature de l’écoulement) 


Illustration de l'intérêt de la méthode : 
Si D, est la traînée mesurée sur une plaque de dimensions Lx h.quand elle est soumise à un 


écoulement de vitesse vı, alors nous avons la relation Eq. (103): 


Di 


AE = Ÿ = = a )oùRe; = Pn 
1 


hı’ R 


L’analyse dimensionnelle par le théorème deBuckinghampermetd’en déduire que pour une 


plaque de dimensions. L, x h: telles que : 


L h 
sl V2 = LA x. Dah l vıhı > Re, = Re, 
2 
x x 


(Similitude de forme)(Facteur d’échelle)(Similitudehydrodynamique) 


et donc IG RE Ge ie 


Dı _ D __ vh _ 
2h2 — 2 p2 D; HE 2h2 D;=D; S D; 
pvini DHU2 TD pvini 








IV-3 Coefficients sans dimension usuels 


En mécanique des fluides, on est souvent amène a manipuler des groupes de variables sans 
dimension, appelés « nombre adimensionnel > ou rapport de similitude >. Ces groupes 
sont construits en faisant des rapports entre des termes apparaissant dans les équations du 


mouvement, ce qui permet de les interpréter physiquement. On distingue ainsi : 


Nombre deReynolds 


L Forces d’inertie 
Re = =>" (104) 


H Forces de viscosité 


avecLune échelle de longueur, vune échelle de vitesse la viscosité du fluide, et psa masse 


volumique. Il exprime le rapport entre les forces d'inertie et les forces de viscositéEq. (104). 
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Il sert notamment à classer le régime d’écoulement en distinguant les écoulements laminaires 
(Re 2000)et les écoulements turbulents (Re 2000). Si on introduit vla viscosité 
cinématique du fluide (v = u /avec pla masse volumique du fluide), alors on a aussi : Re = 


HL, 


Nombre deFroude 


U Forces d'inertie 
Fr = —=—— (105) 
ygL Force de gravité 


% importance pour les écoulements à surface libreEq.(105) 


avecLune échelle de hauteur, v une échelle de vitesse, gl’accélération de la gravite. Le nombrede 
FroudeEq. (105) est le plus souvent interprète comme le rapport de l’énergie cinétique sur l’énergie 
potentielle. Il sert notamment en hydraulique à classer le régime d’écoulement en distinguant 


lesécoulements supercritiques (Fr >1) et les écoulements sub-critiques (Fr <1). 


Remarque : c’est la pesanteur qui est responsable de la forme de lasurface libre : plus Fr est 


grand, moins la surface libre a d’effets surl’écoulement, et inversement. 


Nombre d’Euler 


AP Forces de pression 
Eu = & (106) 

pV Force d'inertie 
Z importance s’il existe de grandes différences de pression au sein de l’écoulement 


Eq.(106). 


Nombre de Mach 
Le nombre de Machexprime le rapport entre la vitesse caractéristique de l'fécoulement et la vitesse 
de propagation du son caractéristique du fluide. Ce nombre quantité donc les effets de 
compressibilitéEq. (107) : 

Ma = U Forces d'inertie 


= -rees d'inertie (197) 


C Force de compression 
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avecuune échelle de vitesse et C la célérité du son (ou célérité des ondes dans l’air). Le 
nombre de Mach sert en aérodynamique ç évaluer la compressibilité de l’air. On distingue 
ainsi les écoulements supersoniques (M >1) et subsoniques (M <1) ; 

st la vitesse du son. 


JPx 


Remarque :C = 


Nombre deStrouhal 


L Forces d'inertielocales 
St = 2 (108) 


V Force d'inertie convectives 


*# importance pour les écoulements non stationnairesEq. (108). 


Nombre de Stokes 


St = + (109) 


avect,le temps de relaxation de la particule (le temps typique de variation de la vitesse quand 
on perturbe l’état d’équilibre de la particule) et le temps caractéristique du fluide (l’échelle de 
temps sur laquelle le fluide s’ajuste a tout changement de la particule). Ce nombreEq. (109). 
sert dans l’étude des écoulements biphasiques (par exemple, une suspension de particules) à 
quantifier les effets biphasiques, c’est-a-dire le couplage entre phases. 

Lorsque St 1, la phase solide est entièrementgouvernée par la phase fluide tandis que pour 


St 1, les deux phases sont découplées. 


Nombre de Prandtl 
Il exprime le rapport entre la viscosité cinématique et la diffusivité thermique du fluide. Il s'agit 
donc d'un nombre caractérisant le fluide. Le nombre de Prandtl Eq. (110). est donc indépendant de 


l'écoulement. 


pratz (110) 
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Nombre de capillarité ou nombre capillaire 


_ Hu 
T 


Ca (111) 


avecuune échelle de vitesse, ula viscosité du fluide, et y la tension de surface. Ce nombre 
Eq.(111)sert à évaluer les effets de tension de surface, par exemple lorsqu’on étale un fluide ou bien 
dans un milieu poreux. Lorsque Ca 1, les effets de tension l’emportent sur les forces visqueuses et 
réciproquement quand Ca 1, la viscosité est tellement grande que les effets de tension de surface a 


la surface libre sont négligeables. 


IV-4 Similitude dans les équations différentielles 


Avant de procéder à l’analyse complète d’un écoulement, il convienttout d’abord de poser les 
hypothèses simplificatrices adéquates.L’évaluation des différents coefficients sans dimension 
relatifs àl’écoulement (Reynolds, Froude...) va en effet permettre de simplifierles équations 


à résoudre. 


Voyons ce qu’il advient de la composante verticale (suivant z) del’équation de Navier- 


StokesEq.(112): 


Introduisons des variables sans dimension : 


u*=u/V x*=x/L p*=p/po 
v*ž=y/V y*=y/L t*=t/t 
w*=w/V ZX) 


oùL, V, po, tsont desgrandeurs caractéristiquesdu système étudié. 
Par extension, on a : 


d _1 à 0? 1 02 
dy Lôx* ôx2 | L? 9x*2 
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d 10 92 1 92 ð 107 


dy  Lðy* a Pay gt rôt 








ô _1 0 32 12? 


ds Ldz* ðz? L ôz 





On obtient alorsl’équation suivanteEq. (113) 


pVôw* pV?( .0w* ðw* ðw* 
L DET 
t ôt* L ð x“ dy* ðz* 

Pu ðP* uV p d2w*  d2w* 


Ron (eee t gps que) 08 (19 


S . pv? 
Divisons toute l’expression®—;nous obtiendronsl’expression suivanteEq. (114) 





L nn: . dw: i TA . dw" 
CES "ox Y” ay " 07 
Pu \ 9P* u \(92w*  d9?2w*  92w* gL 
~ (D) L. ax? T dy"? + z7) US 
avec | j 
wL L AP Po pVL U 
SEE = en 7 = — jp = d 
V Vr pV* pV H JgL 
On peut alors écrire Eq.(115): 
(SE) ðw* pir ðw* pae dw" E ðw* 
ar" ax” ay Y" ar 
- (E PL j ( 1 ) d2w* n d?w* Ñ 9?w* 1 iik 
=- as t e" m ar (Fr) 015) 


Ce qui peut s’interpréter comme suit : 
v si St est très faible : on peut négliger la dérivée instantanée etl écoulement pourra 
être considéré stationnaire. 
v si Eu est très faible : on peut négliger le gradient de pression. 
v si Reest très grand : on peut négliger la viscosité du fluide etl’assimilé à un fluide 
parfait. 


v si Fr est très grand : on peut négliger les effets de la pesanteur 
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IV-5 Application aux maquettes 


L’analyse dimensionnelle, grâce notamment au théorème de Buckingham, permet de 
résumer le comportement d’un système àune relation Eq. (116)entre un nombre restreint de 
grandeurs sans dimension : 


Ti = D(T2, T3, e TEK (116) 


Pour des systèmes complexes la détermination de la fonction ®n ’estaccessible que par 
mesures expérimentales. Ainsi, lors de la miseau point d’un prototype, il est 
économiquement et pratiquement pluspertinent de procéder à ces mesures sur un modèle 
réduit: lamaquette. 

Il faut alors pouvoir transposer les résultats obtenus sur la maquetteà ceux que l’on obtiendra 
sur le prototype : 

On s’arrange donc pour respecter le maximum de similitudes entre lamaquette et le 
prototype. 

A titre d'illustration, reprenons l’exemple de la force de traînéeexercée par un écoulement sur 


une plaque plane Eqs (117, 118): 


Dm Lm Km 
-i = o (Eo => Maquette 117 
nL, " " Lhe "Hafa ha N S 
D L u | 
p p p 
—— = DL K ——— > Prototype (118) 
Dpp (re PpVphp 


On commence alors par se fixer un facteur d’échelle : Y = ASe 
P 
: Lm Lp 
On respecte ensuite le facteur de forme : — = — 
hm bhp 
Si on utilise le même fluide, on a : Pm = Pp €t Um = Up 


Et par conséquent, respecter la similitude de Reynolds revient à : 


h Y 
haa = Rphy Sn He 


Dans ces conditions Eq.(119): 
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Dm D, 


—— = — D - daz ep 5 D, =D 
Da Da hi Do ghg ” SÑ T 


li 


p (119) 


al/a 


Ce résultat, spécifique au problème étudié, montre que par un choixapproprié de la vitesse, la 
mesure expérimentale de la traînée sur unemaquette géométriquement semblable donne 


directement la traînée àlaquelle on doit s’attendre sur le prototype. 


En pratique, tout n’est pas si simple, dans la mesure où l’on ne peutsouvent pas respecter 


simultanément toutes les similitudes. 
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Les unités du système international : Equation aux dimensions, symboles. 


A- Dimension 


Il est commode, pour une grandeur physique X d’introduire sa dimension qui 


seranotée [X]. 


Longueur : L, Masse : M, Temps : T, Intensité : I,  [X]=M°L”T" 


moment d'une force 
pression 

viscosité dynamique 
viscosité cinématique 
tension superficielle 


chaleur, enthalpie 
Entropie 
Conductivité thermique 


Coefficient global 
d'échange thermique 


Capacité thermique 


NB : 





dimensions 
LT! 
LT? 


hertz (Hz) 
newton (N) 


kg/m 


joule (J) 
watt (W) 
Nm 
pascal (Pa) 
Pa.s 

m?/ 
kg s? 


S 


e Les unités du système international SI, unités de base, dimensions, règles d’écriture. 
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B-Alphabet grec 


— Aa | a | 
TS M e 
ES | ë | 
ES | o 





C- Préfixes du SI 


Facteur Préfixe Symbole 
1024 yotta Y 
1021 zetta 
1018 exa 


1015 peta 


10? giga 
10$ méga 


Z 

E 

P 

1012 tera T 
G 

M 

108 kilo k 
h 


102 hecto 
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